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^ ■ Abstract 

X5 ■ 

Qj , In this paper we study the strong gravitational lensing scenario where the lens is a Reissner- 

' Nordstrom black hole. We obtain the basic equations and show that, as in the case of Schwarzschild 

black hole, besides the primary and secondary images, two infinite sets of relativistic images are 
formed. We find analytical expressions for the positions and amplifications of the relativistic images. 
The formalism is applied to the case of a low-mass black hole placed at the galactic halo. 
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O ■ 1 Introduction 
(N ■ 

I The theory of General Relativity predicts the deflection of light in presence of a mass distribution. 
A. Einstein in 1936 t noted that an image due to the deflection of the light of a background star 
I . by another star can have a great magnification if the observer, the lens, and the source are highly 
^Jjl aligned. He also pointed out that the angular separation of the images was too small to be resolved 



X 



by the optical telescopes available at that time. 



It was the discovering of quasars in 1963 which opened the possibility of really observing gravita- 
te ■ tional lensing effects. Quasars are very bright objects, located at cosmological distances, and have a 
central compact optical emitting region. When a galaxy is interposed in the line of sight to them, the 
resulting gravitational magnification can be large and the images are well separated in some particular 
cases. In 1979 the first example of gravitational lensing was discovered (the quasar QSO 0957-1-561 
A,B). 

The weak field theory of gravitational lensing, developed, among others, by Y. G. Klimov, S. 
Liebes, S. Refsdal, R. R. Bourassa, and R. Kantowski, has been successful in explaining the astronom- 
ical observations up to now. This theory is based on a first order expansion of the small deflection 
angle (for a detailed treatment see 0, and references therein). 
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When the lens is a very compact object (e.g. a black hole) the weak field approximation is no 
longer valid. Virbhadra and Ellis recently studied the strong field situation P]. They obtained the lens 
equation using an asymptotically flat background metric and analyzed the lensing by a Schwarzchild 
black hole in the center of the Galaxy using numerical methods. Besides the primary and secondary 
images, they found that there exist two infinite sets of faint relativistic images. Fritelli et al. -C found 
an exact lens equation without any reference to a background metric and compared their results with 
those of Virbhadra and Ellis for the Schwarzchild black hole case. Bozza et al. |H1 obtained analytical 
expressions for the positions and magnification of the relativistic images using the strong field limit 
approximation . 

In this paper we study the lensing situation when the lens is a slowly rotating Kerr-Newman black 
hole. The introduction of charged bodies in the strong-field lensing theory is justified since charged 
black holes are thought to be final result of the catastrophic collapse of very massive (M > 35 Mq) 
magnetized stars. Although selective accretion from the surroundings would neutralize the charged 
black hole if it is located in a high-density medium, there remains the possibility that if the Kerr- 
Newman hole is surrounded by a co-rotating, opposite charged magnetosphere, it might not discharge 
so quickly. Such a configuration would present zero net charge from infinity and consequently it could 
survive for a significant time span (lO^-lO^ yr) if located in a low density environment 6 . Moreover, 
the magnetosphere would have observational effects due to particle acceleration in electrostatic polar 
gaps, similar to those presented by pulsars The relativistic wind created by the Kerr-Newman 
black hole could be responsible for detectable gamma-ray emission [7j. Hence, the study of other 
possible observational signatures from charged black holes presents particular interest. 

In Sec. 2 of this paper we introduce the basic equations using the Reissner-Nordstrom's metric 
for the black hole and a fiat background metric. The Reissner-Nordstrom case can be used also to 
study slowly rotating Kerr-Newman black holes. Highly rotating objects break the spherical symmetry 
introducing unnecessary complications in the lensing calculations, which do not lead to qualitatively 
different results. In Sec. 3 we use the strong field limit approximation to obtain analytical expressions 
for the positions and amplifications of the relativistic images. In Sec. 4 we calculate the positions and 
magnifications of the primary and secondary images using the weak field approximation. In Sec. 5 we 
apply the formalism to a small black hole (7 solar masses) in the galactic halo. We close with Sec. 6, 
where some conclusions are drawn. 



2 Basic equations 

In this paper we use geometrized units (speed of light in vacuum c = 1 and gravitational constant 
G = l). 



Black holes are characterized uniquely by M (mass), Q (charge) and S (intrinsic angular momen- 
tum) [S]. Written in the t, r, 9,(p coordinates of Boyer and Lindquist, the Kerr-Newman geometry has 
the form: 

ds"^ = -^(dt - asin^ Bdipf + ^^[(r^ + a'^)dip - adtf + ^dr^ + p^dO^, (1) 

A 

where 

A = - 2Mr + a"^ + Q^, (2) 

p2 = ^2 _^ ^2 pQg2 0^ 
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Figure 1: Lens diagram. The observer (O), the lens (L), the source (S) and the image (I) positions 
are shown. Do\, Dqs, Dis are, respectively, the observer- lens, the observer-source and the lens-source 
distances, a is the deflection angle and J is the impact parameter. 



The Kerr-Newman geometry is axially symmetric around z axis. The horizon of events is placed 

at 

ru = r+ = M + ^/ M"^ - Q"^ - a? . (5) 

A non-rotating black hole corresponds to an isotropic black hole with charge Q. The metric is the 
Reissner-Nordstrom's one: 

= - (l - ^ + 5) dt" + r2 (sin2 ed^'' + dO^) + (^1 - ^ + ^) ' dr\ (6) 

The horizon is located at 



rn = M + v/m23q2, (7) 

and the photon sphere radius is at 



V = ^M(l + ^lJ|i). (8) 

We study the lens situation shown in Fig. ^ The background space-time is considered asymp- 
totically flat, with the observer and the source immersed in the flat space-time region, which can be 



embedded, if necessary, in a Robertson- Walker expanding Universe. From Fig. 1 we see that the lens 
equation can be expressed |3| 

tan (3 = tan — — ^ (tan 6 + tan(Q — 6)) , (9) 

where /3 and 9 are respectively the angular source and image positions and a is the deflection angle 
due to the black hole. 

Following Sec. 8.5 of 0, we have that the deflection angle for a light ray is 

a{ro) = 2 / , ^ — - vr, (10) 



TO J \ TO Tq J \ ^ I 



where rg is the closest distance of approach. The impact parameter is 

From the lens diagram (Fig. ^) we see that 

J{ro) = Doism9. (12) 
Defining the distances and the charge in terms of the Schwarzschild radius (2M): 

_ j_ _ _ J_ 

^~2M' ^°~2M' 2M' 

, _ -PqI , _ ^os , _ -Pis 

°'"2M' °^"2M' '^"2M' 

= -9- 

^ 2M' 



we have that 



a(xo) = 2 / '^'^ - TT, (13) 



1 q^\-^ 



q 



6(xo)=xo 1 \-^] =do\sin9, (14) 



^H = ^ + y^-g2, (15) 



1 + Jl-^qA. (16) 



The deflection angle a for a light ray passing at the right of the black hole is plotted as a function 
of the closest approach distance xq in Fig. |21 We see that when xq takes values the angle of 
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Figure 2: Deflection angle a (in radians) plotted as a function of the closest approach distance xq. 
The solid line corresponds to Q = 0, the dash-dot line to \Q\ = 0.5M and the dot line to \Q\ = M . 
In each curve the vertical asymptote is placed at xq = Xpg- 



deflection a is greater than 2tt , so the light ray can take several turns around the black hole before 
reaching the observer. In this way, besides the primary and secondary images (with |a| < 27r), we 
have two inflnite sets of relativistic images, one produced by clockwise winding around the black hole 
(a > ) and the other by counter-clockwise winding (a < 0). These images are located, respectively, 
at the same side and at the opposite side of the source. This is qualitatively shown in Fig. |3I 

For a given source position /3 , we must solve Eq. with Eqs. (|13j) and (|14j) to obtain the 
positions of the images. The trascendental equation Q is hard to solve even numerically. 



The magniflcation for a circular symmetric lens is given by 



sin 13 dj3 



sine de 

Differentiating both sides of Eq. ® and with some algebra, we have 

cos 9 



(17) 



dj3 / cos j3 

'de ^ [cose 



where 




/ da 

\de ~ 



de dxQ de ' 



(18) 



(19) 
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Figure 3: The primary(/p) and secondary (Ig) images are formed when the deflection angle is (in 
modulus) smaller than 2tt. The relativistic images (Ir) are obtained when the deflection angle is 
(in modulus) greater than 27r (the flrst with clockwise winding and the flrst with counter-clockwise 
winding are shown). 



with 



da _ /-"^ 

dxo Jxo 



-x^ (xo - 2g2) +xl{x- 2g2) 



-dx, 



/y*3 





1 1 _L 9^ 



(20) 



obtained by differentiating Eq. ((T^ and 



dxQ 



'7^ 



(2x2 - 3X0 + 2q^) 

obtained by differentiating both sides of Eq. (|14|) and doing some extra algebraic tricks. 



(21) 
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3 Strong field limit 



In this section we shall do some approximations. The first one is that when the source and the 
lens are almost aligned we can replace tan 9 hy 6 and tan f3 by /3. For the relativistic images with 
clockwise winding of the rays around the black hole, we can write a = 2mr + Aa^ with n integer and 
< Attn << 1, so that tan(Q — 9) can be approximated by Aon — 9. As in the case of Schwarzschild 
black hole lensing, if a ray of light emitted by the source S is going to reach the observer after turning 
around the black hole, a must be very close to a multiple of 2ir. 
Then the lens equation takes the form j5|: 



As 

Do. 



AOn = 



-Aan, 



and the impact parameter is 



Ai 
2M 



9 = do\9. 



(22) 



(23) 



The relativistic images are formed when the light rays pass very close to the photon sphere. So it 
is convenient to write the closest approach distance xq as 



where < e << 1. Bozza et al 
angle can be approximated by: 



= Xps + e, (24) 
have shown for the Schwarzschild black hole that the deflection 



-2 In 



^ 2 + ^/3 ^ 
18 ^ 



TT. 



(25) 



Prom Fig. we see that the curves representing the deflection angle have similar form for the 
different values of Q. Therefore, we shall also look for a similar approximation 



a 



-Aln(5e)-7r, (26) 
where A and B are to-be-defined positive numbers which will depend only on q = Q/2M. 

A and B are chosen to satisfy 

hm (aexact - "approx.) =0, (27) 

with Oexact given by Eq. l(T^ ^ and Oapprox. is given by Eq. ((^ . Eq. (|?7j) can be written in the form 

Q^exact 



lim 



a. 



approx. 



a 



approx. 



0. 



A necessary (but not sufficient) condition for this limit is: 

C^exact 



lim 

xo^XpB Oapprox. 



1. 



(28) 



(29) 



As a ^ oo when xq 



X 



pS) 



using the L'Hospital's rule in the limit of Eq. ()29() . it is easy to find A: 



A 



lim 

X{) >^3;ps 



-{xq - Xps) 



da, 



exact 



dxc] 



(30) 



^The integral of Eq. 113|l can be calculated in terms of elliptic integrals and it is done in the Appendix, Eq. 1731 . 
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with ^^^^ given by Eq. (pUj) ^. 

To obtain B, we replace Eg .(1261) in Eq. 1)27^ : 



lim [oexact + ^ In S(xo - Xps) + vr] = 0, 



which, using properties of logarithms can be expressed as 

/ "exact + VT 



lim A In 

3^0 '■^'ps 



B{xo - Xps)exp 



so 



then 



lim B{xq - Xps) exp ( '^^^^'^^ ^ ^ ) = l. 



B = lim 



exp 



-(acxact+7r) 

A 



xo^Xps (xo - Xps) 

A and B are given in Table 1 for some values of Q. 



(31) 

(32) 
(33) 

(34) 



Table 1: Numerical values for the coefficients A and B. See text for explanation. 



\Q\ 


O.IM 


0.25M 0.5M 


0.75M 


IM 


A 2.00000 
B 0.207338 


2.00224 
0.207979 


2.01444 2.06586 
0.21147 0.225997 


2.19737 
0.262085 


2.82843 
0.426782 


The ratio ttexact/aapprox. is 


plotted in 


Fig. H for < e < 


0.05. The 


error using 



ttapprox. IS very 

small, less than 1.5 percent for \Q\ < 0.5M and less than 2 percent for 0.5M < \Q\ < M. 

The impact parameter given by Eq. ()14() can be approximated, using a second order Taylor 
expansion in e, by 

b = C + De^, (35) 

with 



3 + V9-32g2 

C= ^ (36) 

4a/6 - 16^2 + 2v^9 - 32g2 



72 - 256g2 



and 



where we have made use of Eqs. (|161 124() . Inverting the Eq. H35|) to obtain £ 



9 - 32g2 + y/9 - 32g2 ^6 - 16^2 + 2v^9-32g2 



(38) 



b-C 



D 

replacing Eqs. and in Eq. we have 



a^-Aln{B\l^^^^^^ \ -TT. (39) 



^Or, in terms of elliptic integrals, by Eq. (|791 of the Appendix. 
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Figure 4: The ratio aexact /oapprox. plotted as a function of e = xq — Xps- The soHd line corresponds to 
Q = 0, the dash-dot line to \Q\ = 0.5M and the dot line to \Q\ = M. 



The position of the n-th relativistic image can be approximated by a first order Taylor expansion 
around a = 2n7r 

~ 6*° - /On Attn, 

with 

^0 = e{a = 2n7r), 

and 



da 



Inverting Eq. (|39() we have that 



1 f D 



then 



and 



1 f D 



exp 



A 



a=2n-K 



■{a + vr) 



— (2n + l)7r 



C 



+ c 



2D 



Pn 



exp 



— (2n + lK 



From Eq. 



Attn 



(40) 
(41) 

(42) 



(43) 

(44) 
(45) 

(46) 
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replacing it in the lens equation, Eq. (|^^ : 



n + 
1 + 



dls On — On 
dos Pn 

dls 

dos Pn 



dls 

dos Pn 



(47) 



Using that pn then a >> 1^ we can neglect the 1 inside the parentheses to obtain 

the approximate positions of the relativistic images: 



= 00 + 



dls 



(48) 



Note that the second term in the last equation is a small correction on 9^. When the source and 
the lens are aligned, (3 = and we obtain the relativistic Einstein rings with angular radius 9^ = 0^. 



The amplification of the n-th image is given by 



Pr 



(3 dl3 

On dOn 



which, using Eq. gives 



Pn 



^0 _|_ dpsPn ^\ dosPn 
dls 



Pdis 



(49) 
(50) 



neglecting the term with f we have 



^ 1 dps ^0 



(51) 



We see that p^ oc p^^r, so the amplification of the relativistic images is very small unless the observer, 
the lens and the source are highly aligned. 

With a similar treatment for the images formed with counter-clockwise winding of light rays around 
the black hole, the positions of the relativistic images are 



On = -Ol + 

"Is 

and the amplifications are given again by Eq. ((^ . but with opposite parity. 



(52) 



To obtain the total magnification of the relativistic images, we must take into account both sets 
of relativistic images, and sum 



PR 



n=l 



(53) 



ra=l 



which, using that the sum of a geometrical series is a" = for \a\ < 1, gives 

n=l 



PR 



1 das 4.D f D 

+ C/ - 



\f3\ disdl^ AB^ \B^l- i _ g-^-M i ' 



(54) 
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4 Primary and secondary images 



We saw in the last section that relativistic images have a very small amplification unless the aligment 
is high. Now we shall obtain the positions and amplifications for the primary and secondary images 
in the aproximation of very small source angle (/?). 



Very small /? implies, from the lens equation, that a is very small too for the primary (P) and 
secondary images (S*), so we can use the weak field aproximation to obtain the positions and amplifi- 
cations. The first step is to calculate the defiection angle a. Expanding the integrand of Eq. H13|) and 
doing the integration we have that 



a 



2 



15^ 
16' 



-vr 



-vrg 



1 



+ 



1 



(55) 



The charge only introduces a small correction in the second order term of the defiection angle. As usual 
in the weak field aproximation, we shall use a first order expansion, taking q 2/xo and 6 « xq ~ c^oi^- 
Introducing q in the lens equation and inverting it to obtain the positions of the images, we have: 



^1^ 



/32 + 4^1 



/52 + 4^ 



-^E + 2/3, 



(56) 
(57) 



where 
by: 



is the Einstein ring angular radius. Using Eq. (|17j) . the amplifications are given 



\0e 



6e J 



+ 4 



1 

+ 2 



I] +4 



(58) 



(59) 



So, the total amplification //„f of the weak field images is 



Dividing Eq. (j5i|) by (|6fl|) we have that 



dos doi 



\d\sdo\) 



where 



AD 



D 



-12-k/A 



/(<?), 



+ C- 



-Q-n/A 



V2AB^ \B^1- e-8^M 1 - g-^^M 



(60) 



(61) 



(62) 



f(q) takes values between 1 x 10^^ and 3 x 10^^, for \q\ G [0,0.5], so the relativistic images are much 
less amplified than the weak field images. To give an example, the prefactor multiplying f{q) is 
{V> /2M)~'^/'^ , if the lens is halfway between the observer and the source, V ~ Dqs/A which, unless the 
black hole is very massive and or near the Earth makes the amplification extremely small. 
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5 An example: black hole in the galactic halo 



In this section we consider as a lens a black hole with mass M = 7Mq , negligible angular momentum 
{S 0) and charge Q placed at the galactic halo (DqI = 4kpc) JHl- The source is a star with radius 
Rs = Rq located in the galactic bulge (Dos = 8kpc). This provides a model for lensing by black hole 
candidates similar to the case discussed in Ref. and is only intended to give some feeling of the 
numerical values implied in the strong field lensing of charged compact objects. 

The positions and magnifications of the relativistic images for a point source are given, respectively, 
by Eqs. (|481 I52|) and (|5H) . The weak field images are formed near the Einstein angular radius 
9e ~ 2.669 mas (milliarc second). The relativistic images have angular radius of the order of the 
relativistic Einstein rings, the outer one is df and the others approach to 6^. In Table 2 the angular 
positions of the horizon, the photon sphere, the outer Einstein ring and the limiting value of the 
Einstein rings are given for different values of Q. 

Table 2: Angular positions of the horizon, the photon sphere, the outer Einstein ring {Of) and the 
limiting value of the Einstein rings (Of^) for different values of Q. 





\Q\ 





0.1 Af 


0.5M 




M 


Oh 


(mas) 


3.455 X 10"« 


3.447 X 10-« 


3.224 X 10" 


-8 


1.727 X 10-8 


Ops 


(mas) 


5.181 X 10"*^ 


5.171 X 10-8 


4.876 X 10" 


-8 


3.455 X 10-8 


ef 


(mas) 


8.987 X 10"*^ 


8.973 X 10-8 


8.595 X 10" 


-8 


6.957 X 10-8 


aE 


(mas) 


8.977 X 10-8 


8.960 X 10-8 


8.581 X 10" 


-8 


6.910 X 10-8 



The source is extended, so we have to integrate over its luminosity profile to obtain the magnifi- 
cation: 

<-) 

where I is the surface intensity distribution of the source and is the magnification corresponding 
to each point of the source. We take the source as uniform, so 

For both the relativistic and the weak field cases, the amplifications for a point source are proportional 
to 1//3: 

To obtain I we use polar coordinates (-R, ip) in the source plane, with i? = in the optical axis, and 
take the source as a disk D{Rc,Rs) of radius Rs centered in Rq- Then 

TiRl ' ^ ' 

and, using that /3 = R/Dos is the angular position of each point of the source, we have 



,/3s) 



djddip 



^al ^Ps ■ ^^^^ 
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Figure 5: Light curves for a lensing event by a black hole with mass M = 7Mq and charge Q placed 
at the galactic halo. The source is a star of radius Rs = Rq placed at the galactic bulge. Upper 
panel, left: total amplification of the primary and secondary images. Upper panel, right: Q = 0, 
total amplification of the relativistic images. Lower panel, left: |Q| = 0.5 M, total amplification of 
the relativistic images. Lower panel, right: \Q\ = M, total amplification of the relativistic images. 
In all plots, the solid line corresponds to /3o = 0, the dash-dot line to (3o = 0.8/3s and the dot line to 
Po = 213,. 

where D(/?c,/3s) is the disk with angular radius /3s = 5.8 x 10~^mas centered in j3c- The last integral 
can be calculated in terms of elliptic functions 0: 



2Sign[l3, - 
+ {ps + Pc)F 



iPs - Pc)E 
TT Ap,p, 

2' {Ps-Pcf 



TT 



4PsPc 



2' iPs-PcY 



+ 
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where 



F 




Jo 



(69) 



'4>o 



( 



1 



A sin' 



2 



(70) 



are elliptic integrals of the the first and second kind respectively, with the arguments </)o and A. 

In Fig. the light curves for a transit event are shown. The relative velocity between the lens 
and the source is v = 300km/s. The angular position of the center of the source is given by /?c = 



Y /3q + y'^^j 1 where /3o is the minimum angular position of the center of the source corresponding 

to t = 0. We see that the transit event takes only a few hours, with great amplification for the 
weak field images and, as expected, a negligible magnification for the relativistic images. However, it 
is interesting to remark that with the presence of charge the amplification of the relativistic images 
increases respect to the Schwarzschild case. For a charge |Q| = 0.5M we have about 10 % more 
amplification than in a similar black hole with no charge. In the case |Q| = IM the total amplification 
of the relativistic images is about 140 % stronger than the corresponding value for a Schwarzschild 
black hole of the same mass. 

6 Conclusions 

In this paper we developed the formalism for charged black hole lensing. We find the positions and 
magnifications of the relativistic images within the strong field limit. The relativistic images are more 
prominent when there is a close alignment between the observer, the lens, and the source. They are, 
though, very faint, even with complete alignment, but much stronger than in the case of Schwarzschild 
black holes for large charges. The angular separation of the relativistic images is beyond the angular 
resolution of current optical technologies. For local black holes (e.g. in the Gould Belt, a star forming 
region at 600 pc), however, if radio or X-ray background sources can be used, the angular separations 
of fias could be resolved with space techniques. 

The NASA Constellation-X (11) mission, to be launched in 2008, is optimized to study the iron 
Ka line features and will determine the black hole mass and spin for a large number of systems. Still, 
Constellation-X will provide an indirect measure of the properties of the region within a few event 
horizon radii. NASA-planned MAXIM mission |T2], a /x-arcsec X-ray imaging mission, will be able 
to take direct X-ray pictures of regions of the size of a black hole event horizon. Both of these space 
mission will have the ability to give us proofs of black hole existence, and possibly, to distinguish even 
among different black hole solutions. The project ARISE (Advanced Radio Interferometry between 
Space and Earth) will use the technique of Space VLBI. The mission, to be launched in 2008, will be 
based on a 25-meter inflatable space radio telescope working between 8 and 86 GHz JJ^. ARISE will 
provide resolution of 15 /iarc-seconds or better, 5-10 times better than that achievable on the ground. 
At frequencies of 43 and 86 GHz, this would imply resolution of light weeks to light months in distant 
quasars and will complement the gamma-ray and X-ray observations. ARISE, then, could also study 
gravitational lenses at resolutions of tens of /xarc-seconds, and as such, could prove important in the 
detection of charged black holes. However, it is not only resolution what is needed in order to detect 
the relativistic images. The main problem will be that they are very much demagnified. They would 
certainly pose a challenge to observations, if ever possible, enhanced by the fact that the relativistic 
images are closed to the non-relativistic ones, which are much more intense. 
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Appendix: Deflection angle in terms of elliptic integrals 

It is convenient for the calculations to express the integrals of Eqs. H13() and ()20|) in terms of standard 

elliptic integrals. 

Eq. ((T^ can be put in the form 

a{xo) = —^ I r—— -TT, (71) 

where ^ 

P{x)=x^ + xl(l- — + %^ {-x'+x-q^), (72) 

V xy 

is a fourth degree polynomial. The four roots of P{x) depend only on xq and q and they are real for 
3^0 > and |Q| < M. Calling these roots ri > r2 > > r4 , the first one is ri = xq and the others 
have long expressions (not given here for lack of space) that can be obtained with standard software. 
Then we can write the exact integral of Eq. (fT^ as 

a(xo) = G(xo)F(0o,A) -TT, (73) 

with 

G{x,) = , ^"l" (74) 
^l-i + %V(^i-^3)(r2-r4) 

and ^ 

F((/>o,A) = / " (l-Asin2</,)-^#, (75) 

JO 

an elliptic integral of the first kind with arguments 



(ri - r3)(r2 - r^)' 
The derivative of the deflection angle can be expressed: 



ho = arcsin x / — ^ , (76) 
V (n - Ti) 

^ ^ in -r4)(r2 -rg) ^^^^ 



dxo dxo OA dxo ocpo dxo 



and can be written as 



^ = \)—^r( ^ f E{(l)o,X) F{^o,X) sin(2,/>o) \ dX 

dxo ''^'^°'^^dxo+'^^'^°^l^ 2(A-1)A 2A +4(A-l)Vl-Asin^0oy ^ 

+ , ^, (79) 

VI - A sin^ (t)Q dxo 
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where 

r<t>o 1 
E{(l)o, A) = / (1 - A sin^ 0) ^ dcj), (80) 







is an elliptic integral of the second kind with the arguments cpo and A given above. 
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